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Abstract 



, We consider the effect of a random magnetic field in the convective zone of the Sun 

' on resonant neutrino spin- flavour oscillations. The expected signals in the different 

experiments (SK,GALLEX-SAGE,H°rnestake) are obtained as a function of the level 
of noise, regular magnetic field and neutrino mixing parameters. Previous results 
obtained for small mixing and ad-hoc regular magnetic profiles are reobtained. We 
■ find that MSW regions are stable up to very large levels of noise (P=0.7-0.8) and they 

, are acceptable from the point of view of antineutrino production. For strong noise 

any parameter region (Am 2 , sin 2 28) is excluded: this model of noisy magnetic field is 
not compatible with particle physics solutions to the SNP. 



1 Introduction 



^h. A neutrino transition magnetic moment can account, both for the observed deficiency 

of the solar neutrino flux and the time variations of the signal. The overall deficit is 
caused by the suppression of the neutrino energy spectrum. The time dependence may 
be caused by time variations of the convective magnetic field of the Sun [U . In transverse 
magnetic fields, neutrinos with transition magnetic moments will experience spin and 
flavor rotation simultaneously (resonant spin-flavor precession, RSFP). The observation 
of electron antineutrinos from the Sun would lead to the conclusion that the neutrinos are 
Majorana particles. There are however stringent bounds on the presence of solar electron 
antineutrinos coming from solar 8 B neutrinos ||[ ||. 

Magnetic fields measured on the surface of the Sun are weaker than those ones within 
the interior of the convective zone. The mean field value over the solar disc is of the 
order of 1 G and in the solar spots magnetic field strength reaches 1 kG. In magnetic 
hydrodynamics, MHD, one can explain such fields in a self-consistent way if these fields 
are generated by the existence of a dynamo effect at the bottom of the convective zone. 
In this region the strength of small scale regular magnetic fields could reach a value of 
100 kG. These fields propagate through the convective zone and photosphere decreasing 

Presenting work made on collaboration with V.B. Semikoz. To be published in the proceedings of 2nd 
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in strength while increasing in the scale giving traces in from of loops in bipolar active 
regions (solar spots). 

Large-scale toroidal magnetic field created by a dynamo mechanism in convective zone 
typically have strengths even less than small-scale r.m.s fields near the bottom of convective 
zone. This is the main reason why one should consider neutrino propagation in the random 
magnetic field of the Sun. The ratio of the r.m.s. random field and the regular (toroidal) 
field may be about ~ 20-50, therefore the problem of RSFP neutrino propagation in noisy 
magnetic field seems to be important. Estimations for the ratio of rms fields to regular 
field are necessarily very roug h. In textbooks @, |, |, § we find the conservative ratio 
r] = (B 2 ) /Bq ~ 1. In more elaborate models the ratio of magnetic energy densities is given 
by the magnetic Reynolds number, (B 2 )/B 2 ~ R^, which may be much bigger than unity 
for plasma with large conductivity. Here a > is a topology index ||, ||, |8[ 

The effect of random magnetic fields in RSFP solutions to the solar neutrino prob- 
lem, SNP, and antineutrino production has been explored previously for simplified models 
[10, 11]. In this work we will deal with the complete problem, we will present calculations 
of neutrino spin flavor conversions in presence of matter and magnetic field. The magnetic 
field will have two well differentiated components. The first component will be a theoreti- 
cally motivated solar magnetic field profile: solution to the static magnetic hydrodynamic 
equations. As a second component, we include a layer of magnetic noise generated at the 
bottom of the convective zone, we have justified that the level of noise in this region can 
certainly be very high. 



2 The Master equation 



We consider conversions u e L,R ^aL,R, a = fi or r, (for definiteness we will refer to fj, 
for the rest of this work) for two neutrino flavors obeying the following master evolution 
equation 



id t 



(VeL \ 
VeR 








S 2 S 





-Ve ~ C 2 5 


-nBX(t) 


S 2 5 


-t*B+(t) 


+ c 2 5 




S 2 5 






fiB+(t) \ 
s 2 5 


-Vv + cidJ 



(VeL \ 
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where c 2 = cos 29, s 2 = sin(2#) and 5 = Am 2 /AE are the neutrino mixing parameters; 
fx = [i\ 2 is the neutrino active-active transition magnetic moment, 



Bf(t) = B± ± (t)+B*(t) 



(2) 



is the magnetic field component which is perpendicular to the neutrino trajectory in the 
Sun. The quantities V e ^{t) are the neutrino vector potentials for u e L and v^l in the Sun 
given by the abundances of the electron and neutron components and by the SSM density 
profile f| . 

The transverse magnetic field B±(t) appearing in Eq. (|2]) is given by B ± (t) = B x (t) ± 
iB y (t) = \B±(t)\Exp(±i$(t)). The nature and magnitude for the regular (Bq) and chaotic 
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parts (B) of the magnetic field will be the subject of the next section. In our calculations 
we will consider always the product pB. Expected values of B pa 1 — 100 kG in the 
Sun convective zone and p = lO -11 /^ would give an expected range for the product 
fiB pa 10~ 8 - 10 -6 p B G ^ 5.6 10~ 17 - 10~ 15 eV or in the practical units which will be 
used throughout this work pB pa 0.1 — 10.0 pi\B^. 

3 Solar magnetic fields 
3.1 Random magnetic fields 

The r.m.s. random component J (B 2 (t)) can be comparable in magnitude with the regular 
one, Bo(t), and maybe even much stronger than Bq, if a large magnetic Reynolds number 
R m leads to the effective dynamo enhancement of small-scale (random) magnetic fields. 

Let us give some simple estimates of the magnetic Reynolds number R m = lv/u m in 
the convective zone for fully ionized hydrogen plasma (T 3> In ~ 13.6 eV ~ 10 5 K). 
Here I ~ 10 8 cm is the size of eddy (of the order of magnitude of a granule size) with 
the turbulent velocity inside of it v ~ va ~ 10 5 cm/s where va = Bo/^/Airp is the Alfven 
velocity for MHD plasma, Bq is a large-scale field in the zone and p is the matter density 
(in g/cm 3 . 

The magnetic diffusion coefficient, v m = c 2 /4Tia con d, appears in the diffusion term of 
the Faraday equation: 

= rot[v x B(t)} + u m AB . (3) 

ot 

Here c is the light speed. The conductivity of the hydrogen plasma is (J con d = ^i/^^ep, 
where u p i = Aire 2 n e /m e = 5.65 x 10 \/n^ s _1 is the plasma (Langmuir) frequency; v ep = 
50n e /T 3 / 2 s _1 is the electron-proton collision frequency, the electron density n e (= n p ) 
(cm -3 ) and the temperature T (K). 

Thus we find that the magnetic diffusion coefficient v m ~ 10 13 (T/1 K)~ 3 / 2 cm 2 s -1 does 
not depend on the charge density n e and it is very small relative to the product v I for hot 
plasma T > 10 5 K » 1 K. From the comparison of the first and second terms in the r.h.s. 
of the Faraday equation Eq. (||) we find that tv/i > f m // 2 , or v m <^ vl ~ 10 13 cm 2 s -1 since 
T/l K 3> 1. This means that the magnetic field in the Sun is mainly frozen-in. Neglecting 
the second term in Eq. (j3|) and using the Maxwell equation rotE = —c~ 1 (dB/dt) we obtain 
the condition for frozen-in field: the Lorentz force vanishes, ~ (E + [vx B]/c) ~ but the 
current j = a con d(E + [v x B]/c) remains finite if the conductivity is large, a C ond — > °°- 

The magnetic Reynolds number 

R m = lvL0 2 pl /(c 2 u ep ) ~ Iv X 10 -13 (T/1 Kfl 2 crr^s" 1 

is huge if we substitute the estimate Iv ~ 10 13 cm 2 s _1 given above. A large value for 
the Reynolds number is a necessary condition for the existence of an effective dynamo 
enhancement in the convective zone. 
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The random magnetic field component in the Sun, B(t), will be described in general 
by an arbitrary, phenomenological, correlator of the form (B(t)B(t )) = (B 2 )f(t — t ). 
We will assume that the strength of the r.m.s. field squared (B 2 ) = i]Bq is parametrized 
by the dimensionless parameter n = > 1 , which it can be, in general, much bigger 
than unity, 77 3> 1. Actually, the estimation of the quantity 77, the ratio of rms fields 
to regular field, for the solar convective zone (and other cosmic dynamos) is the mat- 
ter of current scientific discussions. The most conservative estimate, simply based on 
equipartition, is 77 = constant. According to direct observations of galactic magnetic field 
presumably driven by a dynamo, 77 ~ 1.8 fll3| . A more developed theory of equipartition 
gives 77 ~ 4tt In R m (see Ref. [0]). Let us note that this estimate is considered now as 
very conservative: more detailed theories of MHD turbulence yield estimates like n ~ \/R m 

§!■ 

The correlator function f(t) is unknown a priori but it takes the particular ^-correlator 
form f{t) = LiQ5{t) if the correlation length (for two neighboring magnetic field domains) 
is much less than the neutrino oscillation length, Lq <C lose- Lq can be considered a free 
parameter ranging in the interval 10° — 10 4 km. In the averaged evolution equations it only 
appears the product ijLq. Thus in what follows we will present our results as a function 
of the quantity P which is a simple function of such a product: 

P = i(l + exp(- 7 )), 7 = ^n 2 At^ ^nL (uB ) 2 At. (4) 

The reason for using the quantity P is that it is a good approximation for the depo- 
larization that the presence of noise induces in the averaged neutrino density matrix. At 
is the distance over which the noise is acting. We have supposed in our computations 
that the noise is only effective in a thin layer with thickness At = 0.1 Rq starting at 
r = 0.7 Rq, the bottom of the convective zone. In Table (1) of Ref.(|l5|]) the quantities 

J(B 2 ) and 77 are computed for a given P supposing the reasonable value Lq = 1000 Km. 

For example for P = 0.95, \J (B 2 ) = 10 and n(fiB = 1) = 100 (all uB are given in /J11-B4 
units). 

3.2 Regular large-scale magnetic field in the Sun. The twisting field 

In this work, we will apply for the neutrino conversions described by our master equation 
Eq. (|l]) the self-consistent model of large-scale regular field given in Ref. (jOJ. In this 
model, the global solar magnetic field is the axisymmetric equilibrium solution of the 
MHD static equations (quiet Sun) in the spherically symmetric gravitational field of the 
Sun. The reasonable boundary condition Bq = on the photosphere (r = Rq) is imposed 
in addition. Any field solution to these equations and boundary conditions is a twisting 
field with an, arbitrary, small or large number of revolutions along radius (k = 1,2, • • •, the 
twist rate can be taken as a label to distinguish particular solution within the family). Full 
expressions for the spherical components of the magnetic field can be found in Ref. ( ]i~5[ ) . 
There is only a free parameter in this model. The constant K is related with central field 
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Bcore by: K = B core /2(l — aR & / sin(xR ). The modulus of the perpendicular component 
is of the form: Bq± = B core sin 9f(r)/r where f(r) is some known function of gentle 
behavior. According to this model, the expected magnetic field at the core is typically 
only 2-3 times (or less) the magnetic field at the convective zone. For the values that we 
will consider later, Bq y ~< 100 — 200 kG, the values corresponding at the core are well 
below astrophysical bounds derived from traces of these fields at solar surface. 



4 The averaged master equation 

The master Equation ([!]) can be written in terms of the density matrix p(t) as: 

id t p = [H reg ,p] + pB x (t)[V x , p] + fxB y (t)[V y ,p\. (5) 

The elements of the matrices Ho, V x , V y can be read off the Eq. (g). The B x ,B y are the 
Cartesian transversal components of the chaotic magnetic field. Vacuum mixing terms 
and matter terms corresponding to the SSM density profile given before and the regular 
magnetic part Hamiltonian are all included in H reg . In particular, the matrices V x ^ y are 



Gamma-like matrices in terms of the Pauli matrices a\ t 2 [15|. 

It is our objective in this section to write the differential evolution equation for the 
average density matrix (p). We assume that the components B x ,B y are statistically 
independent, each of them characterized by a delta-correlation function {{B x (t)B y (t')) = 
0). 

(B x>y (t)B x>y (t')) = (Bl y )L 5(t-t). (6) 

In addition, we will make an equipartition assumption for each of the three cartesian 
components B xyz . The averaged evolution equation is a simple generalization (see Ref. 
|p7[ for a complete derivation) of the well known Redfield equation [|l^, [l9|, [2^] for two 
independent sources of noise and reads (Q 2 = Lop 2 (B 2 ) /2 = t]Lq(pBq) 2 /3): 

idt( P ) = [H reg , ( P )] - in 2 [v x , [v x , ( P )]] - in 2 [v y , [v y , ( P )]]. (7) 

It is possible to write the Eq. (0) in a more evolved form. Taking into account the 
particular form of the matrices V x , y and performing a rescaling of the density matrix given 
by: (p(t)) = exp(— AQ 2 t)(p' (t)) we finally obtain the desired averaged evolution equation: 

id t (p') = [H reg , {(/)] + i2Q 2 (V x (p')V x + V y (p')V y ) . (8) 

It is useful however to consider the solution to Eq. (||) when H reg = 0. This is the 
appropriate limit when dealing with extremely low Am 2 or very large energies, for an 
extreme level of noise or when the distance over which the noise is acting is small enough 
to consider the evolution given by H reg negligible. In any other scenario it can give at least 
an idea of the general behavior of the solutions to the full Eq. (|8|). When Hq = only 
the two last terms in the equation remain and an exact simple expression is obtainable by 
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ordinary algebraic methods. The full 4x4 Hamiltonian decouples in 2x2 blocks. If Pf ; i are 
the final and initial probabilities (at the exit and at the entrance of the noise region) their 
averaged counterparts fulfill linear relations among them, schematically: Qf' B = MQf' B 
with Q A,B any of the two dimensional vectors Q A = ((P(u e L — > v e i)), (P(v e L — ► *Vr))), 
Q B = {(P{v e L ^eR))i {P{v e L u hl))) and the Markovian matrix M: 

M = (j P p '/) (9) 

with P defined in Eq. (|j). It can be shown that in this simple case P is exactly the 
final polarization of the density matrix (one of the eigenvalues of the matrix M is equal 
to 2P — 1). In the general case with a finite H reg it can be shown numerically that the 
quantity P still gives a reasonable approximation (< 10%) to the real polarization, at least 
for the cases of interest in this work. 



5 Results and Discussion 



The present status of the Solar neutrino problem in Tables (1) of Refs. [15, 21]. For 
this data, we have calculated the expected neutrino signals in the Homestake, Ga-Ge 
and (Super)-Kamiokande experiments. For this objective, the time averaged survival 
and transition probabilities have been obtained by numerical integration of the ensemble 
averaged master equation (|^) for a certain regular magnetic profile B$± (see Ref. [15] for 



details). The free parameters of our model are four: 5 = Am 2 /2E, s z 2 = sin 29, a 
noise strength parameter (P) and the product of magnetic field and moment (/i-Bo± at 
r = 0.7R Q , 9 S = 7°). 

In Figs. we present the (Am 2 , sin 2 29) exclusion plots from a combined x 2 analysis 
of the three experiments. First we comment the results in complete absence of noise (P=l) 
which are represented in Fig. [I] (Left block of four plots). For negligible or low regular 
magnetic field we observe the high squared mass difference solutions proportioned by the 
matter MSW effect. As the magnitude of the magnetic field increases new solutions appear 
and disappear in a complicated manner. The low angle solution however disappears at 
high magnetic fields, after experiencing some distortion coming from its merging with 
newborn magnetic solutions (compare low angle allowed regions in Figures (b) and (c)). 
The antineutrino production (dashed lines) is in general low and Kamiokande bounds 
are not specially restrictive except at very high magnetic fields. Note in Figure (d) the 
very different behavior of the two existing allowed regions: while the in MSW region 
the antineutrino production is in the 0.1 — 1%, compatible comfortably with Kamiokande 
bounds, the RSFP solution reach a value well above 10% and is excluded by them. It 
seems apparent that there are acceptable particle solutions to the SNP even for very large 
regular magnetic fields. 

The pattern of the electron antineutrino probability is very different when a small level 
of noise (P=0.95, Fig. |l|, right block of four plots) is switched on. For Am? > 10~ 6 eV 2 : 
the antineutrino iso-probability lines follow the characteristic MSW triangular patterns 



6 



Figure 1: The parameter regions consistent with the Homestake, Kamiokande and com- 
bined Gallium experiments in absence of noise (Left Block of four figures): P = 1.00. 
(C.L.= 90%, 95%, 99% (from darkest to lighter shaded areas)). Figures (a,b,c,d) cor- 
respond to values of the regular magnetic field at the bottom of the convective zone: 
\xB = 1,5,8,20 //11-B4. The electron antineutrino averaged probability, is represented by 
the dashed lines. Present Kamiokande bounds are P v jt^ < 0.05. Right Block: the same 
for P = 0.95. 

in this region. The structure of the allowed regions remain unmodified. The electron 
antineutrino yield in these regions is below the 5% level similarly as before. For stronger 
levels of noise (P=0.8, Figs. ||) the same comments can be said. The structure and 
position of the allowed regions from combined total rates are practically unmodified but 
the antineutrino yield impose strong restrictions. For an antineutrino probability smaller 
than 3% only some residual, 90% C.L., allowed regions exist at very small mixing angle, 
Am 2 ~ 10 -6 — 10 -7 eV 2 and moderately high regular magnetic field [Figure (c)]. 

The same regions are still acceptable for P=0.7, (Figs. ^). For this level of noise some- 
thing unexpected happens at extremely high regular magnetic field [Figure (d)] (probably 
too high to be acceptable on astrophysical grounds): a new, large, acceptable region ap- 
pear for Am 2 as 10" 5 . This region disappears again for extremely high chaotic fields 
(P=0.55, Figure not shown here). Note that even for this value of P some residual regions 
with a 90% C.L. are marginally acceptable from reconciliation of all experiment total rates 
and antineutrino bounds (2| if the regular magnetic field is ~ 200 kG (for fi = lO -11 //^). 



Figure 2: Left Block of four Figures: As Fig. g for P = 0.80. Right Block: P = 0.70. 



6 Conclusions 

We have presented calculations of neutrino spin flavor precession in presence a layer of 
magnetic noise at the bottom of the convective zone together with matter and a theoret- 
ically motivated solar magnetic field regular profile. We have justified that the level of 
noise in this region can be very high. We find that MSW regions (Am 2 ~ 10 -5 eV 2 , both 
small and large mixing solutions) are stable up to very large levels of noise (P=0.7-0.8) but 
they are acceptable from the point of view of antineutrino production only for moderate 
levels of noise (P ~ 0.95). The stronger r.m.s field occurs at the convective zone, the 
wider (5m 2 , sin 2 29) region should be excluded when considering the constrain imposed 
by existing antineutrino bounds. For strong noise, P = 0.7 or bigger and reasonable 
regular magnetic field, any parameter region (Am 2 , sin 2 29) is excluded. This model of 
noisy magnetic field is not compatible with particle physics solutions to the SNP. One is 
allowed then to reverse the problem and to put limits on r.m.s field strength, correlation 
length and transition magnetic moments by demanding a solution to the SNP under this 
scenario. 
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